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O I ABSTRACT 

. We present photometric observations of four transits in the WASP- 17 planetary system, ob- 

tained using telescope defocussing techniques and with scatters reaching 0.5 mmag per point. 
Our revised orbital period is 4.0 ± 0.6 s longer than previous measurements, a difference of 
6.6(7, and does not support the published detections of orbital eccentricity in this system. We 
1 model the light curves using the JKTEBOP code and calculate the physical properties of the 

J-j . system by recourse to five sets of theoretical stellar model predictions. The resulting plane- 

tary radius, = 1.932 ± 0.052 ± 0.010 Rj up (statistical and systematic errors respectively), 
provides confirmation that WASP- 17 b is the largest planet currently known. All fourteen 
planets with radii measured to be greater than 1.6Rj up are found around comparatively hot 
(T e s > 5900 K) and massive (Ma > 1.15 M Q ) stars. Chromospheric activity indicators are 
available for eight of these stars, and all imply a low activity level. The planets have small or 
zero orbital eccentricities, so tidal effects struggle to explain their large radii. The observed 
dearth of large planets around small stars may be natural but could also be due to observa- 
tional biases against deep transits, if these are mistakenly labelled as false positives and so not 
followed up. 

Key words: stars: planetary systems — stars: fundamental parameters — stars: individual: 
WASP- 17 



* Based on data collected by MiNDSTEp with the Danish 1 .54 m telescope y Royal Society University Research Fellow 
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1 INTRODUCTION 

Ongoing surveys for transiting extrasolar planets (TEPs) have de- 
tected an unexpectedly diverse set of objects, such as HD 80606 b, 
a massi ve planet on an extr emely eccentric orbit (e = 0.9330 ± 
0.0005: lH6brard et alJ2O10h : super-E arths on very sho r t-period or- 
bits l ike CoRoT-7b and 55Cnce dLeger et al.l 120091 ; IWinn et al.l 
1201 ll) : WASP-18b with a m ass of 10 M Jup and an orbital period of 
0.94 d dHellier et al]|2009l) ; WASP-33b, a very hot planet revolv- 
ing ar ound a metallic-lined pulsating A star dCollier Cameron et al.l 
2010); a system o f six planets transiting the star Kepler- 11 
dLissauer et alj201 l|); and a planet or biting the eclipsing binary star 



system Kepler-16 dDovle et al.ll20Hl) . 

Among these objects, WASP- 17 b stands out as both the 
largest known planet a nd the first found to follow a retrograde orbit 
(Ande rson et alJ[2O10L hereafter A10). However, the reliability of 
its radius measurement was questionable for two reasons. Firstly, it 
rested primarily on a single high-quality transit light curve, whereas 
it is widely appreciated that correlated noise can afflict individ- 
ual light curves whilst remaining un detectable in isolation (e.g. 
iGillon et al]|2007l ; lAdams et ai]|201 ll) . Correlated noise is clearly 
visible in the residuals of the best-fitting model for this transit 
(fig. 2 in A10). Secondly, the orbital eccentricity, e, was poorly con- 
strained, and this uncertainty in the orbital velocity of the planet has 
major implications for the interpretation of the transit light curves. 

The WASP- 17 discovery paper (A10) presented three mea- 
surements of the planetary radius, Rb, based on models with dif- 
ferent assumptions. The preferred model (Case 1) was a straight- 
forward fit to the available data, yielding Rb = l - 74ln' ^ Rj up 
and e = 0.129±°' Jgg. The test of Lucy & Sweenevl dl97l|) . which 
accounts for the fact that a measured eccentricity is a biased esti- 
mator of the true value, indicates a probability of only 83% that 
this eccentricity is significant. Case 2 incorporated a Bayesian 
prior on the stellar mass and radius to encourage them towards 
a solution appropriate for a main-sequence star, and resulted in 
R b = 1.51 ± 0.10 R. Tu p and e = 0.237+0;^. The third and final 
model, Case 3, did not use the main-sequence prior but instead en- 
forced e = 0, and yielded Rb = 1-97 ± 0.10 Rj up . The measured 
size of the planet is clearly very sensitive to the treatment of orbital 
eccentricity. 

iTriaud et alj d20ld) and iBavliss et all d20ld) subsequently 
confirmed the provisional finding that WASP- 17 b has a retrograde 
orbit, from radial velocity observations obtained during transit. 
ITriaud et al.l d2010h also obtained improved spectroscopic param- 
eters for the host star. They assume d e = and found Rb = 
1.986t ; ^ R.T„ P . IWood etal.ld201ll) have detected sodium in the 
atmosphere of WASP-17b, using echelle spectroscopy obtained 
outside and during transit. 

I Anderson et alj d201ll hereafter All) used an alternative 
method to constrain the orbital shape of the WASP- 17 system: mea- 
surements of the time of occultation of the planet by the star in 
infrared light curves obtained by the Spitzer satellite. To first or- 
der, the orbital phase of secondary eclipse (occultation) depe nds on 
the pr oduct e cos to where to is the longitude of periastron dKopall 
Il959h . All obtained ecoso; = 0.00352 ± 0.00075 and e = 
0.02810'ojI, finding ecoso; to be significantly different from zero 
at the 4.8cr level. Inclusion of the Spitzer data, alongside existing 
observations, led to the measurement Rb = 1.991 ± 0.081 Rj up . 
This was achieved without making assumptions about e or cj, 
so represents the first clear demonstration that WASP- 17 b is the 
largest planet with a known radius. One remaining concern was 
that the orbital ephemeris of the system had to be extrapolated to 
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Figure 1. Light curves of WASP-17. The error bars have been scaled to give 
x\ = 1.0 for each night, and in some cases are smaller than the symbol 
sizes. 



the times of the Spitzer observations, potentially compromising the 
measurement of the phase of mid-occultation and therefore e cos uj. 

In this work we present new photometry of three complete 
transits of WASP- 17 b, obtained using telescope-defocussing tech- 
niques. These lead to a refinement of the orbital ephemeris, shed- 
ding new light on the possibility of orbital eccentricity in this sys- 
tem. They also allow a new set of physical properties of the system 
to be obtained, which are more precise and no longer dependent on 
the quality of a single follow-up light curve. We use these data to 
confirm the standing of WASP- 17 b as the largest known planet, at 
R h = 1.932 ±0.053 R Jup . 



2 OBSERVATIONS AND DATA REDUCTION 

We observed four transits of WASP-17 through a Bessell R filter, 
in the 2011 observing season, using the 1.54 m Danish TelescopJ]] 



1 For information on the 1.54 m Danish Telescope and DFOSC see: 

http : //www . eso . org/ sci/ f acilit ies / las ilia/ tele scope s/dlp5 / 
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Table 1. Log of the observations presented in this work. N oha is the number of observations and 'Moon ilium.' is the fractional illumination of the Moon at 
the midpoint of the transit. 



Transit 


Date of 


Start time 


End time 


N obs 


Exposure 


Filter 


Airmass 


Moon 


Aperture 


Scatter 




first obs 


(UT) 


(UT) 




time (s) 






ilium. 


radii (px) 


(mmag) 


1 


2011 04 28 


03:48 


10:10 


128 


120 


R 


1.19 -> 1.00 1.57 


0.212 


32, 46, 65 


0.560 


2 


2011 05 28 


00:58 


03:55 


63 


120 


R 


1.31 -» 1.01 


0.199 


28, 38, 58 


0.762 


3 


2011 06 11 


23:51 


06:51 


152 


100 


R 


1.46 -> 1.00 -> 1.48 


0.826 


30, 42, 60 


0.528 


4 


2011 06 26 


23:32 


05:54 


152 


100 


R 


1.26 -> 1.00 -> 1.45 


0.184 


30, 40, 60 


0.475 



Table 2. Excerpts of the light curve of WASP-17. The full dataset will be 
made available at the CDS. 



BJD(TDB) 


Diff. mag. 


Uncertainty 


55679.66520 


0.000732 


0.000572 


55679.93036 


0.000675 


0.000721 


55709.54791 


0.000352 


0.000761 


55709.67092 


0.019863 


0.000766 


55724.50116 


-0.001133 


0.000560 


55724.79836 


0.000359 


0.000568 


55739.48740 


0.000324 


0.000508 


55739.75282 


0.000740 


0.000499 



semble formed from between two and four comparison stars. Si- 
multaneously to optimisation of the comparison star weights, we 
fitted low-order polynomials to the out-of-transit data in order to 
normalise them to unit flux. A first-order polynomial (straight line) 
was used when possible - such a function is preferable as it is in- 
capable of modifying the shape of the transit - but a second-order 
polynomial was needed for the third transit to cope with slow vari- 
ations induced by the changing airmass. 

The final light curves have scatters in the region 0.5 mmag 
per point, which is close to the best achieved at this (or any other) 
1.5 m telescope dSouthworth et alj|2009d [2010). They are shown 
individually in Fig.Q]and tabulated in Table[2] 



at ESO La Silla, Chile (TableO. Our approach was to defocus the 
telescope and use relativ ely long exposure times of 100-120 s (see 
ISouthworth etaUl2009dlbl) . This technique results in a higher ob- 
serving efficiency, as less time is spent on reading out the CCD, and 
therefore lower Poisson and scintillation noise. It also greatly de- 
creases flat-fielding noise as several thousand pixels are contained 
in each point spread function (PSF), and makes the results insensi- 
tive to any changes in the seeing during an observing sequence. We 
autoguided throughout each sequence in order to keep the PSFs on 
the same CCD pixels, which reduces any remaining susceptibility 
to flat-fielding noise. The second of the observing sequences suf- 
fered from clouds from shortly before the midpoint to after the end 
of the transit, and we were not able to obtain reliable photometry 
from the affected data. 

Several images were taken with the telescopes properly fo- 
cussed, in order to check for faint stars within the defocussed PSF 
of WASP-17. The closest star we detected is 6.9 mag fainter and 
separated by 69 pixels (27") from the position of WASP-17. Light 
from this star did not contaminate any of our observations. 

Data reduction wa s performed as in previous papers 
(Southworth et al. 2009a. b), usin g a pipeline written in IDL0 and 
calling the DAOPHOT package ((stetson 1987) to perform aperture 
photometry with the APEF0 routine. The apertures were placed by 
hand (i.e., mouse-click) and were shifted to follow the positions 
of the PSFs by cross-correlating each image against a reference 
image. We tried a wide range of aperture sizes and retained those 
which gave photometry with the lowest scatter compared to a fitted 
model. In line with previous experience, we find that the shape of 
the light curve is very insensitive to the aperture sizes. 

Differential photometry was obtained against an optimal en- 

2 The acronym IDL stands for Interactive Data Language and is a 
trademark of ITT Visual Information Solutions. For further details see 

|http : / /www . it t vis . com/ProductServices/IDL . aspx| 

3 APER is part of the ASTROLIB subroutine library distributed by NASA. 
For further details see http://idlastro.gsfc.nasa. gov/. 



3 LIGHT CURVE ANALYSIS 

The analysis of our light curves was performed identically to the 
Homogeneous Studies approach established by the first autho r. Full 
details can be found in lSouthwortr] {2008, 2009, 201( 11201 lh . Here 
we restrict ourselves to a summary of the analysis steps. 

The light curves we re modelled using the JKTEBOlfl code 
( ISouthworth et al.ll2004al lbh. The star and planet are represented 
by biaxial spheroids, and their shape governed by the mass ratio. 
We adopted the value 0.0004, but our results are extremely insen- 
sitive to this number. The salient parameters of the model are the 
fractional radii of the star and planet (i.e. absolute radii divided by 
semimajor axis) r\ and rt,, and the orbital inclination, i. The frac- 
tional radii were parameterised by their sum and ratio: 



as the latter are less strongly correlated. 
3.1 Orbital period determination 

Our first step was to obtain a refined orbital ephemeris. Our own 
four datasets were fitted individually and their errorbars rescaled 
to give xt ~ 1-0 with respect to the best-fitting model. This 
step is necessary because the uncertainties from our data reduc- 
tion pipeline (specifically the APER algorithm) tend to be underes- 
timated. We then re-derived the times of mid-transit for the three 
datasets which covered complete transits. Monte Carlo simulations 
were used to assess the uncertainties of the measurements, and 
the resulting errorbars we re doubled based on previous experience 
dSouthworth et all2012allb|) . 

We augmented our three times of mid- transit with 13 measure- 
ments from A10, of which 12 are from observations with Super- 
WASP jPollacco et all 200 6) and one is from their follow-up Euler 

4 JKTEBOP is written in FORTRAN77 and the source code is available at 

http : // www .astro . keele . ac . uk/ jkt /codes / jktebop . html 
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Figure 2. Plot of the residuals of the timings of mid-transit of WASP-17 versus a linear ephemeris. Some error bars are smaller than the symbol sizes. Timings 
obtained from SuperWASP data are plotted using open circles, and other timings (Danish and Euler telescopes) are plotted with filled circles. 



Table 3. Times of minimum light of WASP-17 and their residuals versus 
the ephemeris derived in this work. 
References: (1) A10; (2) This work. 



Time of minimum 


Cycle 


Residual 


Reference 


(HJD - 2400000) 


no. 


(HJD) 




53890.54887 ± 0.00430 


-188.0 


0.01843 




53905.48227 ±0.00380 


-184.0 


0.00989 




53920.42277 ± 0.00250 


-180.0 


0.00846 




53965.23767 ±0.00350 


-168.0 


-0.00246 




54200.57117 ±0.00310 


-105.0 


-0.00449 




54215.52227 ±0.00190 


-101.0 


0.00468 




54271.55737 ±0.00280 


-86.0 


0.00751 




54286.49447 ± 0.00580 


-82.0 


0.00267 




54301.45167 ±0.00570 


-78.0 


0.01793 




54331.32347 ±0.00650 


-70.0 


0.00585 




54555.43697 ±0.00440 


-10.0 


-0.00972 




54566.65087 ± 0.00580 


-7.0 


-0.00227 




54592.80046 ±0.00038 


0.0 


-0.00107 




55679.82838 ±0.00046 


291.0 


0.00084 


2 


55724.65322 ±0.00056 


303.0 


-0.00013 


2 


55739.59522 ±0.00028 


307.0 


-0.00007 


2 



light curve (Table[5J- Taking the time of this follow-up dataset as 
the reference epoch, we find the ephemeris: 

To = BJD(TDB) 2 454 592.80154(50) + 3.7354845(19) x E 

where E represents the cycle count with respect to the reference 
epoch. The reduced x 2 of the fit to the timings is rather large at 
Xu = 2.37, and this has been accounted for in the errorbars above. 
The dominant contribution to this xi arises from the SuperWASP 
timings, which confirms the caveat from A 10 that these may have 
optimistic errorbars. 

A plot of the residuals of the fit (Fig.O at first sight suggests 
the possibility of transit timing variations as might arise from the 
light-time effect induced by a body on a wider orbit around WASP- 
17 Ab, but a periodogram of the residuals from Table[3]does not 
show any peaks above the noise level. We therefore proceeded un- 
der the reasonable assumption that the orbital period is constant. 

3.2 Orbital eccentricity 

As emphasised in Sect.[TJ the possibility of an eccentric orbit is 
an important consideration for WASP-17. The radial velocity mea- 
surements of the star do not strongly constrain eccentricity; the ra- 
dial velocity curve of the star has an amplitude not much greater 
than the size of the errorbars on the individual measurements. 



The observed shape of the transit is not useful because it covers 
only a very small p hase interval (a n essentially ubiquitous situ- 
ation for TEPs; see I Kipping 120080 . The only precise constraint 
on orbital shape was obtained by All, from two occultations ob- 
served using Spitzer. They found the phase of mid-occultation to 
be 0.50224 ± 0.00050, allowing a detection of a non-zero e cos w 
at the 4.8a level. 

Our revised period is 4.0±0.6 s larger than that found by Al 1, 
a difference of 6.6<r, which affects the phase of mid-occultation. 
The actual occultation times are not given by All, but an effec- 
tive time can be inferred from the dates of the observations and 
the orbital ephemerides utilised. We performed this calculation and 
then converted the result back into orbital phase using our new 
ephemeris. This procedure incorporates the necessary conversion 
from the UTC to the TDB timescales. We found the phase of occul- 
tation to be 0.50066, which is consistent with phase 0.5 at about the 
la level. The Spitzer results can no longer be taken as evidence of 
an eccentric orbit in the WASP-17 system. This emphasises the im- 
portance of accompanying occultation measurements with transits, 
in order to avoid uncertainties in propagating ephemerides from 
different observing seasons. 

To confirm this result, we obtained a time measurement which 
represents the actual times of the Spitzer occultations by repeat- 
ing the analysis by Al 1. We found 2454949.5422 ± 0.0016 on the 
HJD(UTC) timescale. After converting to the TDB timescale this 
equates to the phase 0.50059 ± 0.00043, which is equivalent to 
an ecoscj of only 0.00093 ± 0.00068. This differs from zero at 
the 1.4<7 level, which we do not regard as convinc ing evidence of 
orbital eccentricity (see also lAnderson et aT]|2012l) . Further obser- 
vations with Warm Spitzer would be useful in confirming the phase 
of mid-occultation of WASP- 17 b. 



3.3 Light curve modelling 

We modelled the three complete transits from the Danish Telescope 
simultaneously using the JKTEBOP code. The partially-observed 
transit provides a confirmation of the transit depth, but is less reli- 
able than the other ones and has little effect on the solution, so was 
not included in further analysis. The xt of the fit to the three light 
curves is 1.22, which indicates that they don't agree completely on 
the transit shape. Such a situation may arise from astrophysical ef- 
fects such as starspot activity (which can cause changes in the tran- 
sit depth without altering the shape), instrumental effects such as 
correlated noise in the photometry, and analysis effects such as im- 
perfect transit normalisation. Importantly, our possession of three 
independent transits mitigates against all three eventualities, mak- 
ing the resulting solutions much more reliable than those based on a 
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Table 4. Parameters of the fit to the light curves of WASP-17 from the JKTEBOP analysis (top lines). The parameters adopted as final are given in bold. 
Alternative parameters with various constraints on orbital eccentricity and orientation are included, labelled with the e cos lo value adopted. The parameters 
found by other studies are shown in the lowest part of the table. Quantities without quoted uncertainties were not given by those authors, but have been 
calculated from other parameters which were, e cos u and e sin uj values are given to show explicitly the measurements or assumptions relevant to each 
analysis. 



Source 


ecoskL) 


esinoj 




k 






n> 


Danish data 


0.0 assumed 


0.0 assumed 


0.1616 ±0.0021 


0.1255 ±0.0007 


86.71 ±0.30 


0.1436 ±0.0018 


0.01802 ±0.00030 


Euler data 


0.0 assumed 


0.0 assumed 


0.1744 ±0.0080 


0.1322±0.0012 


85.46 ±0.83 


0.1540 ±0.0069 


0.0204 ±0.0011 


Danish data 


0.036 ± 0.033 


-0.10±0.13 


0.1 80 ±0.023 


0.1254 ± 0.0007 


85.9 ±1.1 


0.160 ±0.021 


0.0200 ± 0.0026 


Danish data 


0.00352 ± 0.00075 


-0.027 ±0.019 


0.1576 ± 0.0040 


0.1254 ± 0.0007 


86.87 ± 0.30 


0.1400 ±0.0035 


0.01757 ± 0.00047 


Danish data 


0.00093 ± 0.00068 


-0.027 ±0.019 


0.1591 ±0.0045 


0.1255 ± 0.0007 


86.81 ±0.32 


0.1414 ± 0.0040 


0.01774 ±0.00054 


Adopted solution 


0.0 assumed 


0.0 assumed 


0.1616 ±0.0021 


0.1255 ±0.0007 


86.71 ± 0.30 


0.1436 ± 0.0018 


0.01802 ±0.00030 



A10(Case 1) 


036+ 034 

U.UOD_ Q Q31 


-0.10±0.13 


0.1446 


i293+ a0011 
u.izt)o_ 0014 


87.8ti-° 


0.1281 


0.01658 


A10 (Case 2) 


QS4+ 025 

U.UD<i_ Q 024 


n 9 oo+0.071 
u.zoo_ Q 070 


0.1275 


12Q4+ - 0010 
U,1Z3 -0.0011 


88.161°;™ 


0.1129 


0.01459 


A10 (Case 3) 


0.0 assumed 


0.0 assumed 


0.1622 


l 2 q<5+ - 0010 
u.izao_ Q Q010 


87.8+jo 


0.1436 


0.01855 


Triaud et al. (2010) 


0.0 assumed 


0.0 assumed 


0.1657 


12Q29+ 00077 
u.±z»z»_ 00061 


86.63+°*° 


1467+ 0033 


n ni8Q7+ 00051 
0.01897_ 00040 


All 


00'«2 +0 ' 00076 
u.uuooz_ Q 00073 


-0 027+ - 019 

U - UZ '-0.015 


0.1605 


0.1302 ±0.0010 


86.83!°,;*? 


0.1420 


0.01847 
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Figure 3. Phased light curves of WASP-17 from the Danish telescope (up- 
per) and the Euler telescope (lower), compared to the best fits found using 
JKTEBOP and the quadratic LD law. The residuals of the fits are plotted at 
the bottom of the figure, offset from zero. 

single transit. The excess xt causes larger errorbars to be obtained 
for both error anlysis algorithms (see below) so is accounted for in 
our final results. 

Light curve models were obtained using each of five 
limb darkening (L D) laws, of which four are biparametric (see 
ISouthworthll2008l) . The LD coefficients were either fixed at the- 
oretically predicted valued or included as fitted parameters. We 
found that fitting for one LD coefficient provided a significant im- 

5 Theoretical limb darkening coefficients were obtained by bilinear in- 



provement on fixing both to their theoretical counterparts, but that 
fitting for both led to ill-conditioned models with no further im- 
provement in the quality of fit. We therefore adopted the fits with 
the linear LD coefficient fitted and the nonlinear LD coefficient set 
to its theoretical value but perturbed by ±0.05 on a flat distribution 
during the error analyses (corresp onding to case 'LD-fit/fix' in the 
nomenclature o f lSouthworthll201ll) . This does not cause a signifi- 
cant dependence on stellar t heory because the two LD coefficients 
are very strongly correlated JSouthworth et alj|2007 at) . The results 
for the linear LD law were not used as linear LD is known to be a 
poor representation of reality. 

Errorbars for the fitted parameters were obtained in two 
ways: from 1000 Monte Carlo simulations for each solution, and 
via a residual-permutation algorithm. We found that the residual- 
permutation method returned larger uncertainties for k but not for 
other parameters, indicating that red noise becomes important when 
measuring the transit depth. The final parameter values are the un- 
weighted mean of those from the solutions involving the four two- 
parameter LD laws. Their errorbars are the larger of the Monte- 
Carlo or residual-permutation alternatives, with an extra contribu- 
tion to account for variations between solutions with the different 
LD laws. 

We also modelled the Cousins 7-band light curve from the Eu- 
ler Telescope presented by A 10, in order to provide a direct com- 
parison with our results. The LD-fit/fix option was also the best, 
and correlated noise was found to be important for all photometric 
parameters. A plot of the best fit is shown in Fig.[3]and the corre- 
sponding parameters are tabulated in Table[4] The agreement be- 
tween the Danish and Euler data is poor, especially for k. The Dan- 
ish results should be more reliable as they are based on three transits 
obtained using an equatorially-mounted telescope in excellent pho- 
tometric conditions. The Euler data are more scattered, cover only 
one transit with a gap near midpoint, and were obtained from an 
alt-az telescope so suffer from continual changes in the light path 
over the duration of the observing sequence. We therefore adopted 
the Danish results as final. 

Table|4]also shows a comparison between our photometric pa- 



terpolation in stellar T c g and logg using the JKTLD code available from 

|http : / /www .astro . keele . a c ■ uk/ jkt / cocles/jktld. html 
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Table 5. Derived physical properties of the WASP- 17 system. The upper part of the table contains the individual results from this work; in each case </ b 
3.16 ± 0.20ms" 2 , p A = 0.324 ± 0.012 p and T ' = 1755 ±28K. 



(dEB constraint) 



(Claret models) 



(Y 2 models) 



(Teramo models) 



(VRSS models) 



(DSEP models) 



K b (kms" 1 ) 
M A (M ) 
Ra(Rq') 
loggA (cgs) 
M b (M Jup ) 

^b (Rjup) 
Pb (PJup) 

e 

a (PAT) 
Age (Gyr) 



148.38 ± 
1.272 ± 
1.577 ± 
4.147 ± 
0.473 ± 
1.925 ± 
0.0620 ± 
0.0197 ± 
0.05105 ± 



3.72 

0.096 

0.045 

0.015 

0.035 

0.058 

0.0049 

0.0012 

0.00128 



148.71 ±2.88 
1.280 ±0.076 
1.580 ±0.040 
4.148 ±0.014 

0.475 ± 0.033 
1.929 ±0.052 
0.0619 ± 0.0048 
0.0197 ±0.0012 

0.051 16 ± 0.00099 

2.9+?-* 



148.85 ±1.78 
1.284 ±0.046 
1.582 ±0.029 
4. 148 ±0.011 
0.476 ± 0.029 
1.931 ±0.040 
0.0618 ± 0.0047 
0.0196 ±0.0012 



149.78 ±0.41 
1.308 ±0.011 
1.592 ±0.022 
4.151 ±0.011 
0.482 ± 0.027 
1.943 ±0.033 
0.0614 ± 0.0046 
0.0195 ±0.0011 



148.52 ±1.49 
1.275 ±0.038 
1.578 ±0.027 
4. 148 ±0.011 
0.474 ± 0.028 
1.927 ±0.038 
0.0620 ± 0.0046 
0.0197 ±0.0011 



0.05121 ±0.00061 0.05153 ± 0.00014 0.05 110 ±0.00051 



2.7 



+0.4 



2.1 



+0.0 



2.5 



+0.2 



145.89 ± 
1.209 ± 
1.550 ± 
4.140± 
0.457 ± 
1.893 ± 
0.0631 ± 
0.0200 ± 
0.05019 ± 
3.3 



2.26 
0.056 
0.033 
0.012 
0.029 
0.043 
0.0048 
0.0012 
0.00078 

0.6 

o.s 



Table 6. Final physical properties of the WASP-17 system (with statistical and systematic errorbars) compared to published measurements. Eccentricity is 
included to illustrate the difference approaches taken to obtain each set of results. 





This work (final) 


A10(Case 1) 


A10 (Case 2) 


A10 (Case 3) 


Triaud et al. (2010) 


All 


e 


0.0 adopted 


12g +0.i06 


n 9 o 7 +0.068 
U - ZO '-0.069 


0.0 adopted 


0.0 adopted 


O.oa8±g;gi| 


M A (M ) 


1.286 ± 0.076 ±0.020 


1.20 ±0.12 


1.16 ±0.12 


1.25 ±0.13 


1.20±0.12 


1.306 ±0.026 


fl A (R©) 


1.583 ± 0.040 ±0.008 


1 no + 0.20 

l-3»_ 01g 


1 200 +a081 
1,z ""-0.080 


1.566 ± 0.073 


1 57 n+0.067 


1.572 ±0.056 


log g A (cgs) 


4.149± 0.014 ±0.002 


4.23 ±0.12 


4.341 ± 0.068 


4 143+0.032 




4.161 ±0.026 


Pa (pq) 


0.324 ±0.012 


45+ - 23 


0.67+° 3 16 


n o 9 o+0.035 


304+ - 016 
"■ JU, -0 .020 


0.336 ±0.030 


M h (M Jup ) 


0.477 ± 0.033 ±0.005 


4QO+ 059 

U.19U_ 05 g 


"■™-0.060 


4QS+ - 059 
V.1DH_ 05g 


453+0-043 
u.ioo_ 035 


0.486 ± 0.032 


(Rjup) 


1.932 ± 0.052 ±0.010 


I 74+0.26 
1 ' ' -0.23 


1.51 ±0.10 


1.97 ±0.10 


1.986_ 074 


1.991 ±0.081 


9b (ms~ 2 ) 


3.16 ±0.20 


3.63+J* 


5 0+ 1 ' 1 
°- U -0.9 


9 Q9+0-36 
z ' az -0.33 




2.81 ±0.27 


Pb (PJup) 


0.0618 ± 0.0048 ±0.0003 


092+ 054 


144+OO 42 
u - iy ^-0.031 


0.0648_ Q 0090 




0.0616 ±0.0080 


T e ' q (K) 


1755 ± 28 




1557 ± 55 


1756± 2 




1771 ±35 





0.0196 ± 0.0012 ±0.0001 












a(AU) 


0.05125 ± 0.00099 ± 0.00027 


0501 + 0017 
u.uoui_ 001g 


0494+ 0017 
u-u*»*_ 001g 


0507+ 0017 
u - uou ' -0.0018 


0.0500 ±0.0017 


0.05150 ±0.00034 


Age (Gyr) 


9 7 +0.6 +0.6 
z -' -1.0 -0.6 


o Q+0.9 


i 9+2.8 
A,z -1.2 


3 1+ 1 - 1 




2.65 ± 0.25 



rameters and those published by other researchers. The values from 
iTriaud et al.1 feOld) and All are in good agreement with our re- 
sults, except for the parameter k. This is as expected, because those 
authors had only the Euler light curve with which to constrain the 
transit shape. The parameters from A10 are more discrepant, pri- 
marily because their favoured solution is for a large orbital eccen- 
tricity. 

To demonstrate the influence of eccentricity, we have rerun the 
above analyses on the Danish data for three ec centric orbits, ap- 
plying constraints using the method described by Sout hworth et al.l 
j2009ch . For the first set of constraints we used e cos u = 0.036 ± 
0.033 and esincj = -0.10 ± 0.13 (A10), for the second we 
adopted ecoso; = 0.00352 ± 0.00075 and esincj = -0.027 ± 
0.019 (All), and for the third we specified e cos ui = 0.00093 ± 
0.00068 and e sin w = -0.027 ± 0.019 (Sect.HH. 



The results for these three alternative sets of constraints are 
shown in Table|4] It can be seen that k is unaffected but, as ex- 
pected, rA±»"b and thus r\ and r b are very dependent on the treat- 
ment of eccentricity. That large eccentricities can be ruled out is 
vital for precisely measuring the properties of WASP-17. The very 
small e cos ui value allowed by the transit and occultation timings 
(Sect. 13.2^ results in photometric parameters which are close to the 
la errorbars of the zero-eccentricity result, and illustrates the small 
change in the measured system properties to be expected if the 1.4a 
measurement of e cos uj turns out to be real. 



4 THE PHYSICAL PROPERTIES OF WASP-17 

The physical properties of the WASP-17 system can be obtained 
from the adopted photometric parameters (Tabled, t he orbital ve- 
locity amplitude of the star (Ka = 52.7±2.9ms _1 ; Tria ud et al.l 
l201fj) . its effective temperature (T e g) and metallicity ( [§] ), and a 
constraint from theoretical stellar evolutionary models. Full details 
of our approach can be found in I Southw orth (2009, 120101) . 

One immediate problem faced here is the diversity of the 
published T cff measurements of WASP- 17 A: A10 find 655 ± 
100 K from analysis of echelle spectra; pTriaud et aT i 201o|) ob- 
tain 6 650 ± 80 K from similar observations, and iMaxted et al.l 
d201 lh deduce 6500 ± 75 K f rom the InfraRed Flux M ethod (IRFM; 
iBlackwell & Shallij 1 1971 iBlackwell etal] Il980h . We adopted 
?cff = 6550 ± 100 K as this encompasses all three determinations 
but leans towards the IRFM value, which should be the method 
with the least dependence on stellar theory and analysis technique. 
The corresponding metallicity is [§1 = -0.25 ± 0.09 (A10). 

Our approach was to guess a value of the orbital velocity of the 
planet (K b ) and combine it with the photometric parameters and 
Ka to calculate th e physical prop erties of both bodies using stan- 
dard formulae (e.g. lHil ditch 2001). The predicted properties of the 
star for the calculated mass were then found via interpolation in a 
set of theoretical model predictions. K b was iteratively adjusted to 
find the best agreement between the known va and calculated 
and between the measured and predicted T ff values. This was done 
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Table 7. Compilation of selected physical properties of TEP systems containing a planet larger than 1.6 Rj up . The projected spin-orbit misalignment mea- 
surable from the Rossiter-McLaughlin effect is denoted by A. All asymmetric errorbars have been averaged, and the statistical and systematic errorbars have 

been added in q uadrature, in order to fit in to the Table. 

Refere nces : mlHartman et alj feOllh: (2)lHartman et al] J2012b; OllLatham et alj feplOl); (4)ISouthworthl <201 A; (5) lFortnev etail 201 ll); (6) lKonacki etal 
20051) : ffl lSouthworthl l20ld) : (8) lKonacki et alj 120031) - (9) ISnerien et alj 120091) : aOl lMandushev et alj j2007n : (ll) lHebbet al] 120091) : < 12) |joshi et al 



20091) : (13) lEnoch etaU <20 1 lb . 



.System 


1VI A 


«A 


J eff 




6 


• 

Period 




R b 


rpl 

cq 


A t ) 


Inn- f?' 

log H HK 


KclS 






I.M-0J 


(IV 






W 


(M Jup ) (R Jup ) 


I. IV 








HAT-P-32 


1.160 


1.219 


6207 


-0.04 


0.0 


2.150 


0.860 


1.789 


1786 






l 




±0.041 


±0.016 


±88 


±0.08 






±0.164 


±0.025 


±26 








HAT-P-33 


1.375 


1.637 


6446 


±0.07 


0.0 


3.475 


0.762 


1.686 


1782 






l 




±0.040 


±0.034 


±88 


±0.08 






±0.101 


±0.045 


±28 








HAT-P-40 


1.512 


2.206 


6080 


+0.22 


0.0 


4.457 


0.615 


1.730 


1770 




-5.12 


2 




±0.073 


±0.061 


±100 


±0.10 






±0.038 


±0.062 


±33 








HAT-P-41 


1.418 


1.683 


6390 


+0.21 


0.0 


2.694 


0.800 


1.685 


1941 




-5.04 


2 




±0.047 


±0.047 


±100 


±0.10 






±0.102 


±0.064 


±38 








Kepler-7 


1.41 


2.028 


5933 


+0.11 


0.0 


4.885 


0.453 


1.649 


1619 




-5.099 


3.4 




±0.10 


±0.038 


±50 


±0.05 






±0.068 


±0.038 


±15 








Kepler- 12 


1.16 


1.490 


5947 


+0.07 


0.0 


4.438 


0.430 


1.695 


1485 






5.4 




±0.12 


±0.050 


±100 


±0.04 






±0.049 


±0.030 


±25 








OGLE-TR-10 


1.277 


1.520 


6075 


+0.28 


0.0 


3.101 


0.68 


1.706 


1702 






6.7 




±0.083 


±0.100 


±86 


±0.10 






±0.15 


±0.054 


±54 








OGLE-TR-56 


1.34 


1.737 


6119 


+0.25 


0.0 


1.212 


1.41 


1.734 


2482 






8,4 




±0.10 


±0.045 


±62 


±0.08 






±0.18 


±0.059 


±30 








OGLE-TR-L9 


1.43 


1.499 


6933 


-0.05 


0.0 


2.486 


4.4 


1.633 


2034 






9,4 




±0.10 


±0.043 


±58 


±0.20 






±1.5 


±0.046 


±22 








TrES-4 


1.339 


1.834 


6200 


+0.14 


0.0 


3.554 


0.897 


1.735 


1805 


+6.3 ± 4.7 


-5.104 


10,4 




±0.086 


±0.087 


±75 


±0.09 






±0.075 


±0.072 


±40 








WASP-12 


1.38 


1.619 


6300 


+0.30 


n nis+ 024 


1.091 


1.43 


1.825 


2523 




-5.500 


11,4 




±0.19 


±0.079 


±100 


±0.10 




±0.14 


±0.094 


±45 








WASP- 14 


1.35 


1.666 


6475 


+0.0 


0.088 


2.244 


7.90 


1.633 


2090 


-33.1 ±7.4 


-4.923 


12,4 




±0.12 


±0.087 


±100 


±0.2 


±0.003 




±0.47 


±0.092 


±59 








WASP-17 


1.286 


1.583 


6550 


-0.25 


0.0 


3.735 


0.477 


1.932 


1755 


-148.5±f 4 


-5.331 


This work 




±0.079 


±0.041 


±100 


±0.09 






±0.033 


±0.053 


±28 






WASP-48 


1.19 


1.75 


6000 


-0.12 


0.0 


2.144 


0.98 


1.67 


2030 






13 




±0.05 


±0.09 


±150 


±0.12 






±0.09 


±0.10 


±70 









for ages ranging from zero to the point at which the star evolves to 
log g < 3.5, leading to a final set of best-fitting physical properties 
and age for the system. 

Statistical errors in th e input parameters wer e propagated us- 
ing a perturbation analysis ( Southwort h et alj2005l) . Systematic er- 
rors arising from the use of theoretical predictions were estimated 
by running separate solutions for each of five i ndependent stel- 



lar m odel tabul ations: Claret l Claretl 2004]) . Y 2 jPemarque et alj 
2004), Teramo iPietrinferni et al.ll2004l)."T^55 fVandenBer g et al.1 



2006) and DSEP footter et alj2008l) . Finally, a model-independent 



set of results was generated using an empirical calibration of stel- 
lar properties found from well-studied eclipsing binaries, a process 
we label "dEB constrai nt". The empirical calibration follows the 
approach introduced by [Enoch et aL (2010) bu t with the improved 
calibration coefficients derived bv lSouthworthl q201 lh . 

The results for each approach are given in Table|5] The mass, 
radius, surface gravity and density of the star are denoted by Ma, 
Ra, log (?a and pa, and of the planet by Mb, Rb, <7b and p b - The 
orbital semimajor axis is a, T c ' q is the equilibrium temperature of 
the planet (negle cting albedo and heat redistribution) and 6 is the 
Safronovl dl972h number. We find that the agreement between the 
model sets is excellent excep t for the DSEP models, which are quite 
discrepant. ISouthwortblfeOllI) noticed that the agreement between 
models deteriorates at lower metallicities, as is experienced here. 
Our final physical properties for the WASP-17 system (Table[6} 



were therefore calculated from the results obtained using the other 
four model sets (Claret, Y 2 , Teramo and VRSS). Table|6]also con- 
tains results from published studies of WASP-17, which are in good 
agreement overall but contain some optimistic errorbars. 



5 WHAT CAUSES SUCH LARGE PLANET RADII? 

We have measured the radius of WASP- 17 b to be R h = 1.932 ± 
0.053 Rjup (adding the statistical and systematic errors in quadra- 
ture), confirming its status as the largest planet currently known. 
We also find no significant evidence for orbital eccentricity, remov- 
ing an additional contribution to the uncertainties of the measured 
properties of the system. 

Its closest competitor, HAT-P-32 b, suffers a similar problem 
concer ning the influence of orbital shape on the resulting planetary 
radius (Har tman et al.|[2oTTh . However, in its case the possibilities 
are switched because the eccentric-orbit alternative has a positive 
e sin uj compared to the previously-postulated negative e sin u> for 
WASP-17. The preferred solution for HAT-P-32, with a circular or- 
bit, results in a smaller radius (R b = 1.789 ± 0.025 Rj up ) com- 
pared to the eccentric-orbit alternative (R b = 2.037 ± 0.099 Rj up 
withe = 0.163 ±0.061). 

The other twelve TEPs with radii above 1.6Rj up are listed 
in Table|7j with a summary of the physical properties of them and 
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Figure 4. Plot of orbital period versus surface gravity for the known transit- 
ing planets. Planets with radii above 1.6Rj up are shown using black filled 
circles, and WASP-17 b is highlighted with an extra circle. The other planets 
are shown with lighter open circles. 
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Figure 5. Plot of stellar T c r versus planet radius for the known transiting 
systems. Other comments are the same as in Fig|4] 



their parent stars. This list has been boosted by the addition of two 

objects, OGLE-TR-56 and WASP- 14, based on major revisions to 

I 1 J k 

their radii by Southworth (2012 ). The newly-d iscovered systems 

WASP-78 and WASP-79 Ismallev et all 1201 2h are not included 

here as the radii of their planets are uncertain. We now discuss the 

properties of the fourteen planets with radii above 1.6 Rj up . 

Firstly, Fig.|4] shows that their orbital period distributior0 
is not exceptional, and that they are consistent with the known 
correla tion between period and surface gravity ( Southworth et al. 
2007b). The masses of all but two of them (OGLE-TR-L9 and 
WASP-12) are in the interval 0.4-1.0 M Jup . These objects also do 
not represent a high-eccentricity population: al l have orbits which 
are (or are almost) c irculafl Tidal heatin g IBoden heimer et all 
l200ll : ljackson et alj20oll : llbgui & Burrowsi 2009) is therefore not a 
viable proposition to explain their large radii. Orbital misalignmen t 
may be relevant: WASP - 14 is misaligned Jjohnson et al.l l2009h . 
WASP-17 is retrograde jTriaud et al.ll2010h . and TrES-4 is axi- 
ally aligned jNarita ■ et alfeoioh . Thus two of the three planets with 
obliquity measurements are misaligned. 

However, a noticable feature of the fourteen large planets is 
that they orbit host stars with T eff > 5900 K and M A > 1.15 M Q . 
This establishes a connection between a bloated planet and a hot 
star. Fig.|5]and Fig.|6]show that the large planets are preferentially 
associated with hotter and more massive host stars. The associa- 
tion does not work the other way: such stars also possess TEPs 
with smaller radii representative of the general planet population. 
A correlation with host star [2^1 was suspected but not found. An 
important factor in bloating planets above their expected size there- 



6 Data taken from the Transiting Extrasolar Planet Catalogue (TEPCat, 
available at: http : / / www . astro. keele.ac.uk/jkt/tepcat/ 

7 It has been noticed tha t orbital eccentricity is correlated with planet mass 
iSouthworth et al.l20 09c) but plots of planet radius versus eccentricity (not 
shown here) indicate that there is - if anything - a deficit of large-radius 
eccentric planets. 



fore seems to be the T e g of their host star. This may be due to 
the enhanced UV flux from such stars, but such a possibility does 
not explain why less inflated planets are found around stars with 
T cS > 5900 K. 

A similar situation occurs with irradiation: the large planets 
have high T c ' q s (or equivalently large specific incident stellar fluxes) 
of 1600 K or more, but so do many other smaller planets. The con- 
nection between inflated radii and hi gh T e ' q is well-known (see 
Baraffe et ld1 l20ld : lEnoch et alj|20lj and references therein) but 
the simultaneo us existence of small planets with high is not 
yet understood. Laug hiin et al.l d201 ll) quantified the correlation be- 
tween the radius anomaly (observed radiu s versus those predicted 
by the models of IBodenheimer et alj2003f ) of TEPs and their T e ' q s, 
in the form of a power law. Whilst their fig. 2 exhibits apprecia- 
ble evidence for this claim, the large planets discussed here remain 
outliers even in that diagram. 

We have sea rched for values for the chromospheric activity in- 
dicator log R^k l Noves et alii 19841) for the host stars of our large 
planets. lKnutson et alj fedlCh give values of -5.331 for WASP-17, 
-5.099 for Kepler-7, -5.104 for TrES-4, -5.500 for WASP-12 
and —4.923 for WASP-14. These logi? HK values are representa- 
tive of inactive stars, suggesting a correlation between inflated plan- 
etary radii and low chromospheric activity. HAT-P-32 and HAT - 
P-33 possess values of the related S index ( Vaugha n et alJ[l978T) 
which indicate a low activity level , despite their compar atively high 
rotation rates and velocity jitter (Har tman et al.l 120111) . HAT-P-40 
and H AT-P-41 have similarl y quiet log i? H K values of —5.12 and 
-5.04 faartman et ail2012l) . 

Therefore eight of the fourteen stars in question have mea- 
sured activity indicators, and all suggest low chromospheric activ- 
ity. Perhaps the increased high-energy phot on flux from more ac- 
tive stars acts against large planetary radii. iKnutson et alj ^20 id) 
found a correlation with the atmospheric properties of planets, 
in that inactive stars possess planets with temperature inversions 
where as plane t s arou nd active stars do not have inversions. How- 
ever, iHartmanl j2oToh found no correlation between log -Rhk an d 
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Figure 6. Plot of stellar mass versus effective temperature for the known 
transiting planets (data taken from TEPCat). Other comments are the same 
as in Fig[4] The errorbars of the planets have been surpressed for clarity, 
with the exception of those fourteen with a large radius. 
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Figure 7. Plot of stellar versus planet radius for the known TEP systems. 
Other comments are the same as in Fig|4] Points without errorbars repre- 
sent those discovered via the space telescopes CoRoT and Kepler, whose 
discovery biases are very different from those of ground-based surveys. 



planet radius, so the low values for the fourteen host stars in ques- 
tion may be an artefact of their T c g distribution. 

Another possibility to explain the large planetary radii is that 
the more massive stars, around which the larger planets are found, 
have shorter main-sequence lifetimes. They will therefore be on 
average younger than the less massive TEP host stars. The large 
planets could simply be at an earlier stage of their evolution. This 
is, however, at odds with the low activity levels of the host stars, 
which implies that t hey are not particularly young. The models by 
iFortnev et alj d2007l their fig. 5) show that it is possible for planets 
to be 2.0 Rjup or above if they are young (of order 10 7 yr) and low- 
mass (less than 1 Mj up ). These criteria are not satisfied by any of 
our large planets. Explaining the radii of the large TEPs is therefore 
only viable if there is a large systematic error in our estimation 
of the age of their host stars, which is unlikely but certainly not 
impossible. 



5.1 Can observational biases explain the properties of large 
planets? 

An important question is whether the correlation between large 
planets and hot and massive host stars is real, or is it merely a 
phantom arising from observational selection effects? The corre- 
lation could easily be suspected to arise from detection biases as 
a function of transit depth. Fig.|7] shows a plot of the radii of the 
host stars versus their planets. We have overlaid dotted lines to in- 
dicate loci of approximately constant transit depth, calculated using 
JKTEBOP and assuming quadratic LD with both coefficients equal 
to 0.30. The subset of planets discovered via space telescopes are 
down-played in Fig.|7] as these should be much less biased against 
finding systems with shallow transits. It is immediately clear that 
the large planets do not stand out in this diagram as having unusu- 
ally deep, and therefore easily detectable, transits. 

Fig.|7] shows that most known TEPs have transit depths of 
about 1-3%. Firstly, we are straightforwardly biased against find- 



ing small planets around hotter (and therefore bigger) stars, as the 
transit depths in these systems are small. The relative paucity of 
small planets around big stars is plausibly explained by lower de- 
tection probabilities for transits less than 0.5% deep. Secondly, the 
absence of large planets around small stars might be explicable by 
either natural rarity or a bias against deep transits. Such objects 
may have a low follow-up priority within planet-search consortia 
if it is believed that deep transits are associated with false positives 
such as eclipsing binary star systems. As an example, if WASP- 17 b 
orbited an unevolved 0.9 Mg star then the transit depth would be 
roughly 5%. 

None of the observables discussed above pre-determines the 
existence of big planets. They are found only around hotter stars, 
but some such stars possess small planets. They have quite high 
T e ' q s, but so do many smaller planets. They do not give rise to un- 
usually deep or shallow transits, which rules out the more simple 
observational biases, and the metallicities of their host stars are not 
exceptional. The question of what causes their bloated radii remains 
unsettled. 



6 SUMMARY AND CONCLUSIONS 

Whilst WASP- 17 was widely regarded to be the largest known 
planet, its radius was uncertain as it was based primarily on one 
follow-up transit light curve which shows moderate correlated 
noise. In this work we present three new high-precision transit light 
curves, obtained using telescope-defocussing techniques, and im- 
prove the measurement of its radius. 

We have refined the orbital ephemeris of the system using our 
new data, which increase the temporal baseline by 3.1 years. Our 
revised orbital period is 4.0 ± 0.6 s longer than previous measure- 
ments, a difference of 6.6a, and this change is sufficient to bring 
the observed time of occultation (All) into line with that expected 
for a circular orbit. Further observations would allow this result to 
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be checked. In the case of WASP-17, circularity of the orbit favours 
a larger value for the planetary radius. 

We have modelled our new data, along with the single follow- 
up light curve presented in the discovery paper (A10), using the 
JKTEBOP code. We paid careful attention to the treatment of limb 
darkening and to obtaining reliable uncertainties. The physical 
properties of the system were then derived using our new photo- 
metric and published spectroscopic results. Remaining uncertain- 
ties and discrepancies in the existing T e g and \yf\ measurements 
of the host star impose a bottleneck on the quality of the resulting 
properties, and new radial velocity observations would also be use- 
ful in refining the measurement of the planet's mass and therefore 
its density. 

WASP-17 b is the largest known planet, with a radius of i?b = 
1.932 ± 0.053 R.j U p. Another eleven planets are known with radii 
greater than 1.6Rj up . They are found only around comparatively 
hot (Ted > 5900 K) and massive (Ma > 1.15 M ) stars, and have 
correspondingly high equilibrium temperatures (T e ' q > 1600 K 
with the exception of Kepler-12) and equivalently incident fluxes. 
But other stars of similar mass and T e g possess smaller planets, 
whilst other planets with similar T c ' q (or equivalently specific inci- 
dent flux) do not have such enlarged radii. One possible discrimi- 
nating feature is that all eight of the fourteen host stars with mea- 
sured activity indicators are chromospherically inactive. 

The set of fourteen large planets do not have unusual transit 
depths. However, planets of this size around cooler stars may have 
an anomalously low discovery rate if their deep transits (of order 
5%) discount them from detailed follow-up observations. High- 
precision radial velocity measurements are expensive in terms 
of telescope time, so dubious TEP candidates with unexpectedly 
deep transits may be prematurely rejected as false positives. A re- 
evaluation of such objects will either yield scientifically valuable 
discoveries, or dismiss the existence of large planets around small 
stars. 

The fourteen planets with radii greater than 1.6 Rj up all have 
circular (or nearly circular) orbits, so their large size cannot easily 
be attributed to tidal heating. Of the three published measurements 
of the axial alignments of these planets, one reveals a retrograde 
orbit, one a misaligned orbit, and one indicates alignment. Obser- 
vations of the Rossiter-McLaughlin effect on the remaining eleven 
systems could either verify or discount the possibility that axial 
alignment is a relevant aspect of large planets. 
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